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Abstract. We report on magnetic field measurements made 
in the innermost coma of 67P/Churyumov-Gerasimenko 
in its low-activity state. Quasi-coherent, large-amplitude 
(6B/B^1), compressional magnetic field oscillations at 
^40 mHz dominate the immediate plasma environment of 
the nucleus. This differs from previously studied cometary 
interaction regions where waves at the cometary ion gyro- 
frequencies are the main feature. Thus classical pickup-ion- 
driven instabilities are unable to explain the observations. We 
propose a cross-field current instability associated with new¬ 
born cometary ion currents as a possible source mechanism. 

Keywords. Magnetospheric physics (plasma waves and in¬ 
stabilities; solar wind interactions with unmagnetized bod¬ 
ies) - solar physics, astrophysics, and astronomy (magnetic 
fields) 


1 Introduction 


Typically ionization of atoms and molecules of cometary 
origin is the most important process for the interaction 
of strongly outgassing comets and the solar wind. Dur¬ 
ing encounters of the ICE, Sakigake, and Giotto space¬ 
craft with active comets 21P/Giacobini-Zinner, IP/Halley, 
and 26P/Grigg-Skjellerup, large-amplitude plasma waves 
and turbulence have been one of the most pronounced ob¬ 
servational findings in the cometary magnetosphere ([Tsuru- 


|tani and Smith] |1986[ ?; [Neubauer et al.| |1986[ [Glassmei^ 


let al.[[T989l [Glassmeier and Neubau^|1993||Volwerk et al.| 

2014|l. Ion ring-beam instabilities (Wu and Davidson 1972|l 
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and non-gyrotropic phase space density-driven instabilities 
( |Motschmann and Glassmel^ |1993| l are the source mecha¬ 
nism of these waves. In the spacecraft (s/c) frame of refer¬ 
ence (and pickup ion frame), those waves were detected at 
the cometary H 20 ''' ion gyro-frequency. 

Rosetta’s journey ( Glassmeier et al[| 2007a| l alongside 
comet 67P/Churyumov-Gerasimenko now allows for electro¬ 
magnetic waves to be studied at the beginning of cometary 
activity, at the birth of the cometary magnetosphere ( |Nilsson| 
|et al.|[20T5| l. It should be noted that, under these low-activity 
conditions, typical solar wind-cometary interaction regions 
like the bow shock and magnetic pileup region are not ex¬ 
pected (e.g., Koenders et al. 201 3[ Rubin et al. 2014| l and 
also not observed (Nilsson et al. 2015| l. We shall report on 
wave observations at distances of 2.7-3.6 AU from the Sun 
and 10-1000 km from the comet. 


2 Mission and instrumentation 


RPC-Mag_P 67P/C-G CSEQ September 10, 2014 



Figure 1. Example of magnetic field observations made onboard 
the Rosetta spacecraft on 10 September 2014, 23:45-23:55 UTC. 
The position vector of the spacecraft in the comet-centered so¬ 
lar equatorial (CSEQ; for details see text) coordinate system was 
(3.9,-20.6,20.4) km. 


Rosetta arrived at 67P/Churyumov-Gerasimenko on 6 Au¬ 
gust 2014 at a heliocentric distance of 3.6 AU. The spacecraft 
was initially put into a ~ 100 km orbit around the comet’s nu¬ 
cleus. Observations reported here are limited to the dayside 
inner coma and sampled over a time span of a ^4 months 
(from August to November 2014). Using measurements of 
the Rosetta Orbiter Spectrometer for Ion and Neutral Anal¬ 


ysis Cometary Pressure Sensor (ROSINA COPS) (Balsiger 


et al. 2007|l, the cometary activity at this heliocentric dis¬ 


tance was determined to be below 4 x 10^® s“^. This produc¬ 
tion rate is 2-3 orders of magnitude lower than at any other 
previous cometary encounter where pickup ion waves were 
detected (8 x 10^® s“^ at IP/Halley, down to 7 x 10^^ s“^ at 
26P/Grigg-Skjellerup) (see Richter et al. 2011| l. 

The Rosetta orbiter is equipped with a suite of plasma in¬ 
struments - the Rosetta Plasma Consortium (RPC) set of par¬ 
ticle and field sensors dCarr et al.| |2007l l. RPC-MAG, the 
tri-axial fluxgate magnetometer system ([Glassmeier et al. 


|2007b^ , consists of two sensors mounted on a 1.5 m boom, 
separated by 0.15 m. The short boom length implies that 
the spacecraft is heavily contaminating the magnetic field 
measurements. At this stage of the investigation it was not 
possible to completely remove these quasi-static spacecraft 
bias fields from the measured magnetic field values. The dy¬ 
namic range of RPC-MAG is ±16 000nT, and its resolu¬ 
tion 0.03 nT. Although the magnetometer is capable of ac¬ 
quiring the magnetic field measurements with sampling rates 
up to 20 Hz, the data presented here correspond to the in¬ 
strument’s normal operational mode, i.e., a sampling rate of 
1 Hz, which is sufficient for the purposes of our investiga¬ 
tion. The magnetic field observations are represented in a 
comet-centered solar equatorial (CSEQ) coordinate system. 
The +x axis points from the comet to the Sun, the +z axis 
is the component of the Sun’s north pole of date orthogonal 
to the +x axis, and the +y axis completes the right-handed 


reference frame. The origin of the coordinate system is the 
comet’s center of mass. 


The Rosetta Ion and Electron Sensor (RPC-IES) ( jBurchj 
jet al.j 2006| l and the Ion Composition Analyser (RPC- 
ICA) ( [Nilsson et al.||20()6| ) provide information on cometary 
ions produced in the coma of 67P/Churyumov-Gerasimenko. 
Observations of the neutral gas number density made by 
ROSINA COPS are used as well. 


3 Observations 

Upon arrival at the comet on 6 August 2014, RPC-MAG 
started to detect large-amplitude, quasi-coherent magnetic 
field fluctuations. During the aforementioned observational 
period, we were able to collect ^ 3000 cases of wave activity 
with typical frequencies of ~ 40 mHz. Eigure 1 shows an ex¬ 
ample of these waves. The wave activity is clearly visible in 
all three components. Peak-to-peak amplitudes are of the or¬ 
der of 4 nT, which is about twice as large as the ambient solar 
wind magnetic field at this heliocentric distance. The oscilla¬ 
tions are neither purely transverse nor purely compressional. 

A preliminary minimum variance analysis between Au¬ 
gust and November 2014 does not show any preferred di¬ 
rection of wave propagation, with respect to neither the so¬ 
lar wind flow nor the mean magnetic held direction. A full 
discussion of the minimum variance analysis results will be 
presented elsewhere. As only single point observations are 
available, because only the ROSETTA orbiter magnetometer 
was operating and the lander was still attached to the orbiter, 
no information about a typical wavelength can be inferred at 
this time. 

The quasi-coherent nature of these fluctuations is also 
clearly visible in power spectral density distributions 
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67P/CG 2014, OctoberlO 10:40-10:57 UT 



Figure 2. Example of a spectrum of magnetic field fluctuations 
observed in the innermost coma of 67P/Churyumov-Gerasimenko 
on 10 October 2014, 10:40-10:57 UTC (CSEQ position vector 
(—0.3, —10.1, —1.7) km). The power spectrum has been calculated 
by application of a standard fast Eourier transform routine to a 
17 min time interval. A box car window is used to average in the 
frequency domain. Eighteen degrees of freedom are used. The in¬ 
set shows the distribution of peak spectral frequencies in the period 
August-November 2014. The confidence interval is shown in the 
bottom right comer. 


(Fig. ig. The steep spectral slope at frequencies beyond 
the peak frequency is not uncommon for the observations 
at 67P/Churyumov-Gerasimenko; typical spectral slopes are 
between —3 and —5. The wave spectra typically exhibit a 
single-peak structure as shown in Fig.[^ 

The peak frequencies exhibit a Rayleigh-type distribution 
grouped around 40 mHz (Fig.|^. Further statistical analyses 
do not show any clear correlation of the peak frequencies 
with the ambient magnetic field magnitude. Such a correla¬ 
tion would be expected if the observed frequency coincides 
with the local proton gyro-frequency /p (/p ~ 40 mHz with 
B = 2.5 nT). Though the agreement between the peak fre¬ 
quency of the Rayleigh-type distribution and the local gyro- 
frequency is striking, we argue here that the observed waves 
are not in proton cyclotron resonance as a clear correlation 
between magnetic field magnitude and frequency is missing. 


Wave activity was first observed by the RPC-MAG instru¬ 
ment at a day side distance of approximately 100 km from 
the nucleus and steadily increased up to a distance of 30 km 
(Fig.|3j. Hourly trace spectral densities have been integrated 
over the frequency range 30-80 mHz and divided by twice 
the vacuum permeability to determine the magnetic energy 
density of the observed fluctuations. Closer to the nucleus the 
magnetic energy density seems to saturate. The spread of the 
distribution at around 20 and 30km is caused by Rosetta’s 
trajectory with respect to the comet. During the first months 
after arriving at the comet the spacecraft was often been po¬ 
sitioned in the so-called bound orbits with a fixed distance 


to the nucleus, mostly in the terminator plane (Hassig et al. 


|2015| l. Hence, more observations are made at these distances. 
The wave activity variations at constant distance are caused 
by temporal variations in the neutral gas density produced 
by localized gas sources on the nucleus’ surface and mod¬ 
ulated by the comet’s rotation (Hassig et al. 2015| ), leading 
through the ionization of cometary neutrals to temporal vari¬ 
ations in the cometary plasma density. The variations also re¬ 
flect a dependence on solar wind variations as well as eleva¬ 
tion and azimuth angle of Rosetta’s position. Future work is 
planned to separate these dependencies. For the range of 30- 
100 km the spatial variation in the magnetic energy density 
is found to be TF cx with a = 8.04 ± 0.27 (for a 95 % 
confidence interval ol G [7.77; 8.31], yielding a relative error 
of 100 • ^ = 3.3 % ), equivalent to a quartic decrease in the 
wave amplitude with distance. 

Comparing wave intensity with the neutral gas number 
density as measured by ROSINA COPS exhibits a clear 
global relation between both quantities. COPS also detected 
a neutral gas number density, N, above its noise level at 
a distance of about 100 km (Fig. |^. The neutral density 
decreases as W cx with /?= 1.43 ±0.07 (for a 95 % 
confidence interval /? G [1.36; 1.50]). The deviation 100- 
=28.5% from the theoretically expected value 
/3theo = 2 ( [Coates and Jones] [2009 | l is due to Rosetta first 
approaching from the afternoon side (45° phase angle) be¬ 
fore moving to the terminator plane (90° phase angle) and 
therefore encountering different insolation conditions at the 
sub-spacecraft location. 

Here only a global relation is discussed. More detailed 
analyses on the relation between neutral gas density and 
magnetic field variations are currently being prepared; how¬ 
ever, such studies require consideration of the rotation of the 
nucleus with a period of about 12.4 h as well as any inhomo¬ 
geneity in the active regions. 


4 Discussion and possible wave source mechanism 

Comparing the radial variations in both the magnetic energy 
density, W, and the neutral gas density suggests an approx¬ 
imate global relation W cx with 7 = 5.62 ± 0.45 (for a 
95% confidence interval 7 G [5.17; 6.07]). A more detailed 
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Figure 3. Radial variation in the magnetic energy density of the 
low-frequency wave activity in the spectral range 30-80 mHz. The 
inset shows the radial variation in the density of cometary neutrals 
as measured hy the neutral gas monitor COPS of the ROSINA in¬ 
strument. 


correlation analysis between the individual values of both 
quantities at the same radial distance reveals a linear Pear¬ 
son correlation coefficient of r = 0.52. This correlation re¬ 
flects the large spread of both wave energy and neutral gas 
density. As waves are generated locally under time-varying 
solar wind conditions and propagate within the inner coma, 
this weak correlation is not surprising. Because the produc¬ 
tion rate of heavy cometary ions is proportional to the neutral 
gas density, we conclude that wave activity in general is con¬ 
trolled by the cometary ion production rate. 

The observed frequencies concentrate between 25 and 
75 mHz (Fig. |^. In the coma the local H20’*" ion gyro- 
frequency is in the range of 0.8-3.2 mHz. Thus, there are 
clear differences between the observed frequencies at the 
birth stage of the cometary magnetosphere and the heavy 
ion cyclotron frequencies, as expected for well-developed 
cometary interaction regions. This points towards a new gen¬ 
eration mechanism for the type of wave reported here for a 
weakly outgassing comet. 

The size of the innermost interaction region is much less 
than the Larmor radius of the newborn ions. The pickup ion’s 
initial velocity is about 0.8 km s“^; acceleration by the inter¬ 
planetary electric field, E, up to velocities, | v |, of a few tens 
of kilometers per second is observed (Nilsson et al. 2015| l by 
RPC-IES and RPC-ICA. After ionization the newborn ions 
are moving transverse to the ambient magnetic field, B, and 
the solar wind flow in the direction of the electric field, con¬ 
stituting a cross-field electric current density. The cometary 
ion motion is controlled by the electric field. Lorentz forces 
are not yet important in this innermost coma region. The ra¬ 
tio of the electric force to the Lorentz force is of the order of 
the ratio of the gyro-period Tq to the lifetime t of the new¬ 
born ions: \ E\ /\ vxB\Ki Tq/t. Newborn cometary ions 
with r ^ Tg are essentially unmagnetized. Lorentz forces 
become important only on scales comparable to and larger 


than a cometary ion Larmor radius. Within distances to the 
nucleus smaller than a Larmor radius, i.e., in the Larmor 
sphere ( |Sauer et arj|1998| l, physical processes different from 
those generating the classical pickup-ion-related ion cyclo¬ 
tron waves are important. 

A newborn water ion flux, n-v, of at least 10^° m“ 


has been observed by the RPC-ICA sensor (Nilsson et al. 


|2015| l. Due to the sensor’s limited field of view, the actual 
flux will be higher. Assuming a flux density 3xl0i°m-2s-i 
and singly charged ions the electric current density is esti¬ 
mated to be about j ^ 4.8 x 10“® Am“^. Electromagnetic 
instabilities attributable to such a cross-field current have 
been studied in the past (Chang et al. 1990[ Sauer et al. 


|1998| l, but not for conditions typical for the plasma situation 
at 67P/Churyumov-Gerasimenko. Assuming that the cross- 
field current is driven unstable, a transverse wave number 
k± can be estimated using Ampeere’s law: k± = ^qSJ/SB. 
With Sj ^ j and SB « 1 nT, a value of A:_l « 6 x 10“® m“^ 
results, corresponding to a transverse wave length of 524 km. 
It should be noted that this wave number component is trans¬ 
verse to the mean magnetic held as well as the cross-field cur¬ 
rent. The polarization of the observed waves is neither purely 
compressional nor purely transverse. This points towards off- 
angle propagation. We therefore assume that all three com¬ 
ponents of the wave vector are of comparable magnitude. 

The wavelength reported here is larger than the scale of 
the dayside inner coma where wave activity has been ob¬ 
served. However, the generation region is probably signifi¬ 
cantly larger. But as waves generated trough cometary inter¬ 
action at larger distances are buried in the pre-existing solar 
wind turbulence, the signal was not detected by RPC-MAG. 
We need to await a further increase in the cometary activity to 
see a further expansion of that region where cometary waves 
clearly stand out from the solar wind turbulence background. 

The frequency of the unstable mode may be estimated us¬ 
ing the beam-mode dispersion relation ( |Chang et al.| |1990| 
1998[l w « /c|| • w , where fc|| is the beam parallel 


Sauer et al. 


wave vector component and v the cometary ion beam veloc¬ 
ity- _ _ 

Using an observational estimate (Nilsson et al. 20151 of 
the ion velocity ri = 40kms“t, we obtain an angular fre¬ 
quency omega oj ~ 0.24 rads“t, i.e., a wave frequency / ^ 
38 mHz. 

This value is comparable to the observed frequencies. 

It should be noted that the suggested wave source is not 
co-moving with the solar wind flow. The cross-field current 
source is due to freshly ionized cometary ions which are not 
yet moving with the solar wind velocity as they did not have 
time to be accelerated to solar wind speed on the time and 
distance scales we are looking at. The wave source is almost 
fixed in the nucleus frame of reference. Therefore, Doppler 
effects can be neglected as the Rosetta spacecraft is only 
slowly moving with respect to the nucleus (Vg/c ~ 1 m s“^). 

As the wave number is proportional to the current den¬ 
sity perturbation, k (x Sj (x nv, the dispersion relation pro- 
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vides the following approximate expression between wave 
frequency, ion density and ion velocity: / cx nv^. Assuming 
a constant electric field accelerating the particles, the ion ve¬ 
locity should increase linearly with distance. For an ion den¬ 
sity decreasing with the square of the distance, the frequen¬ 
cies should therefore not exhibit any major dependence on 
distance, which is what we observe. Of course, this is only a 
first conjecture stimulating future analysis. 

Our model furthermore allows for the apparent saturation 
of wave activity at around 30 km distance to be understood: 
the newborn ions need to be accelerated to constitute a sig¬ 
nificant current and the waves need to grow. However, further 
detailed theoretical modeling is required to validate the con¬ 
jecture of a cross-field current-driven instability causing the 
newly detected low-frequency wave activity in the Larmor 
sphere of 67P/Churyumov-Gerasimenko. 
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